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Abstract. We show that the presence of a temporal electromagnetic field on 
cosmological scales generates an effective cosmological constant which can account 
for the accelerated expansion of the universe. Primordial electromagnetic quantum 
fluctuations produced during electroweak scale inflation could naturally explain the 
presence of this field and also the measured value of the dark energy density. The 
behavior of the electromagnetic field on cosmological scales is found to differ from 
the well studied short-distance behavior and, in fact, the presence of a non-vanishing 
cosmological constant could be signalling the breakdown of gauge invariance on 
cosmological scales. The theory is compatible with all the local gravity tests, and 
is free from classical or quantum instabilities. Thus we see that, not only the true 
nature of dark energy can be established without resorting to new physics, but also 
the value of the cosmological constant finds a natural explanation in the context of 
standard inflationary cosmology. This mechanism could be discriminated from a true 
cosmological constant by upcoming observations of CMB anisotropics and large scale 
structure. 
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1. Introduction 

The nature of dark energy, which is believed to be responsible for the present phase 
of accelerated expansion of the universe [H [21 El H] , still remains unknown. Despite its 
phenomenological success, the simplest description in terms of a cosmological constant 
(ACDM model) suffers from an important naturalness problem, since the measured value 
of A, corresponding to pa ~ pm ~ (2 x 10^^ eV)^, finds no natural explanation in the 
context of known physics. Moreover, the fact that today matter and dark energy have 
comparable contributions to the energy density, turns out to be difficult to understand if 
dark energy is a true cosmological constant. Thus, the energy density of a cosmological 
constant remains constant throughout the history of the universe, whereas those of the 
rest of components (matter or radiation) grow as we go back in time. Then the question 
arises as to whether it is a coincidence (or not) that they have comparable values today 
when they have differed by many orders of magnitude in the past. Notice also that if A 
is a fundamental constant of nature, its scale (around 10~^ eV) is more than 30 orders 
of magnitude smaller than the natural scale of gravitation, G = Mp'^ with Mp ~ 10^^ 
GeV. On the other hand, if A is just an effective parametrization of dark energy, still 
a proper understanding of the underlying physics would be needed in order to explain 
the measured value. 

Alternative models have been proposed in which dark energy is a dynamical 
component rather than a cosmological constant. Such models are usually based on 
new physics, either in the form of new cosmological fields or modifications of Einstein's 
gravity O El [71 El [9]. However, they are generically plagued by classical or quantum 
instabilities, fine tuning problems or inconsistencies with local gravity constraints. 

In this paper we explore the possibility of understanding dark energy from 
the standard electromagnetic field, without the need of introducing new physics 
(previous works on models of dark energy based on vector fields can be found in 
[ini[IIl[l2l[l3l[ll[l5l|16l[I7]). We will show that the behavior of electromagnetic 
fields on very large (super-Hubble) scales differs from the well studied short-distance 
(sub-Hubble) behavior. Thus, on super-Hubble scales, the time component of the 
electromagnetic field grows linearly in time in the matter and radiation eras, giving 
rise to a cosmological constant contribution in the electromagnetic energy-momentum 
tensor (the potential gravitational effects of longitudinal electromagnetic fields were 
considered in a different context in [18j). At late times this contribution becomes 
dominant giving rise to the accelerated phase. As a possible generating mechanism, we 
calculate the spectrum of super-Hubble electromagnetic modes produced during inflation 
from quantum fluctuations and flnd that the correct value of the dark energy density 
can be naturally obtained in the case in which inflation took place at the electroweak 
scale. 
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2. Cosmological electromagnetic fields 

We start by writing the electromagnetic action including a gauge-fixing term in the 
presence of gravity: 



S = I d'^Xy/^ 



(1) 



167rG 4 2 
The gauge-fixing term is required in order to define a consistent quantum theory for the 
electromagnetic field [19], and we will see that it plays a fundamental role on large scales. 
Still this action preserves a residual gauge symmetry ^ + 9^0 with 00 = 0. 
Einstein's and electromagnetic equations derived from this action can be written 

as: 

- \Rg,. = 8nG {T,, + T^;) (2) 
y^F^"" + AV'^V^A'^ = (3) 

where T^^, is the energy-momentum tensor for matter and radiation and is the 
energy-momentum tensor of the electromagnetic field. Notice that since we will be using 
the covariant Gupta-Bleuler formalism, we do not a priori impose the Lorentz condition. 
The effect of the high conductivity of the universe in the matter and radiation eras will 
be discussed below. 

We shall first focus on the simplest case of a homogeneous electromagnetic field 
(zero mode) in a fiat Robertson- Walker background, whose metric is given by: 

ds^ = dt^ - a{t)'^6ijdx'dx^ (4) 

In this space-time, equations read: 

Ao + 3HAo + 3HAo = (5) 

A + HA =0 (6) 

with H = d/a the Hubble parameter. 

Notice that ([5]) implies that the gauge-fixing term exactly behaves as a cosmological 
constant throughout the history of the universe, irrespective of the background 
evolution. Indeed, for homogeneous fields we have: 

j^{V,An = j^{A, + 3HAo)=0 (7) 

We can solve and (jS]) during the radiation and matter dominated epochs when 
the Hubble parameter is given by H = p/t with p = 1/2 for radiation and p = 2/3 for 
matter. In such the solutions for (El) are: 

Aoit) = A^t + A^t-^P (8) 
A{t) =A+t^-P + A' (9) 

where and are constants of integration. Hence, the growing mode of the temporal 
component does not depend on the epoch being always proportional to the cosmic time 
t, whereas the growing mode of the spatial component evolves as t^/^ during radiation 
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and as t^^'^ during matter, i.e. at late times the temporal component will dominate over 
the spatial ones. 

On the other hand, the (° o) component of Einstein's equations adopts the following 
form: 

SttG 



3 



(10) 



^ Po, + pAo + Pa 

la=R,M 

where R, M stands for radiation and matter respectively and: 

fiAn = A (j^^'-lj + iHAA + lAlj (11) 

= (12) 

Notice that we need A > in order to have positive energy density for Aq. 
Besides, when inserting the solutions ([H]) and ([9]) into these expressions we obtain 
that = Pao^ Pa ~ P'x^"'^ V^A'^ = const as commented before. Thus, the 
field behaves as a cosmological constant throughout the evolution of the universe 
since its temporal component gives rise to a constant energy density whereas the 
energy density corresponding to A always decays as radiation. Moreover, this fact 
prevents the generation of a non-negligible anisotropy which could spoil the highly 
isotropic CMB radiation (see [20] for a more general discussion). Finally, when the 
universe is dominated by the cosmological constant arising from the gauge-fixing term, 
both the Hubble parameter and Aq become constant leading therefore to a future de 
Sitter universe. Let us emphasize that according to (JTj), paq always contributes as 
a cosmological constant. As the observed fraction of energy density associated to a 
cosmological constant today is ^ 0.7, we obtain that the field value today must be 
Ao{to) ~ 0.3 Mp. 

The effects of the high electric conductivity cr can be introduced using the magneto- 
hydrodynamical approximation and including on the r.h.s. of Maxwell's equations the 
corresponding current term, which is given by J^ — J^u^u^ = aF^yU^ with the velocity 
associated to the comoving observers. Notice that the strict neutrality of the plasma, 
which is consistent with a vanishing electric field, implies J^m^ = 0, and finally, the 
current can be written as: = (0, cr{doAi — diAo)). Notice that electric neutrality also 
implies that conductivity does not affect the evolution of v4o(t). The infinite conductivity 
limit simply eliminates the growing mode of A{t) in iQ. The inhomogeneous case, 
corresponding to 7^ modes, will be discussed in next section. 

We still need to understand which are the appropriate initial conditions leading 
to the present value of Aq. In order to avoid the cosmic coincidence problem, such 
initial conditions should have been set in a natural way in the early universe. In a very 
interesting work [21] , it was suggested that the present value of the dark energy density 
could be related to physics at the electroweak scale since pA ~ (^Ivk/-^^')^; where 
Mew ~ 10'^ GeV. This relation offers a hint on the possible mechanism generating 
the initial amplitude of the electromagnetic fluctuations. Indeed, we see that if 



Cosmological electromagnetic fields and dark energy 



5 



such amplitude is set by the size of the Hubble horizon at the electroweak era, i.e. 
AQitEwY ~ ^Ew^ then the correct scale for the dark energy density is obtained. 
Thus, using the Friedmann equation, we find H'^-^^ ~ M^-^^/Mp, but according to 
(fTTj) . ~ H'^Al ~ const., so that paq ~ H^^r ~ {M^^r/Mp)^ as commented before. 

A possible implementation of this mechanism can take place during inflation. Notice 
that the typical scale of the dispersion of quantum field fluctuations on super-Hubble 
scales generated in an inflationary period is precisely set by the almost constant Hubble 
parameter during such period Hj, i.e. (A^) ~ H] [22j. The correct dark energy density 
can then be naturally obtained if initial conditions for the electromagnetic fluctuations 
are set during an inflationary epoch at the scale Mj ~ Mew- Let us make these 
arguments more precise. 



3. Quantum fluctuations during inflation 



We shall look at the electromagnetic perturbations generated during inflation in order 
to determine its primordial power spectrum. In this case it is more convenient to 
use conformal time rj defined by means of dt = adrj and to introduce the conformal 
components of the field = {aAo, A). Besides we shall focus on a single Fourier mode 
of the vector field with wave vector k and decompose the field in temporal, transverse 
and longitudinal components with respect to k. In this frame, equations ([3]) read: 



A" - 



2n' + 



Ok 



2ik 



1 + A 
2A 



A[. - HAw 







(13) 



A" 



A'.' 



_!_;. + /C A±k 

k^XAu 



aaA', 



±k 



2ikX 



1 + A 
2A 



Ok 



aa{A'uf, — ikAok) 



with ' = and Ti. = aH is the Hubble parameter in conformal time. He have included 
for completeness the current term on the right-hand side as commented before. Notice 
that once again the electric neutrality of the universe implies that the evolution equation 
for the temporal component is not modified. During inflation the electric conductivity 
of the universe is negligible and this term can be safely neglected so that in the following 
we shall set a = 0. 

It is easy to see from equations ( fT3l) that the transverse modes are just plane waves 
irrespective of the expansion rate. On the other hand, the components Aok and A\\k are 
coupled to each other even in the absence of gravity. This is due to the fact that we are 
working with arbitrary A and not using the simple Feynman gauge A = — 1. 

Let us first consider quantization in Minkowski space-time, with Ti = Ti' = 0. The 
decomposition in Fourier modes can be written as follows: 



^0 



d^k 



+ 



2A;o(27r)3 
.1 + A 
^1^ 



i^—^(sio + a||)A;?7 + ao ) e ^''^ 



1 - A 
aj" + Sit)kr] + a^ 



ikx 
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d^k 



2ko{2 



1 + ^. 

t- — ^(ao + a||)fc?7 + a| 



-ikx 



+ 



1 — A 



Akx 



(14) 



with fco = 1^1 = k. Now, in order to have the canonical commutation rules 
[A^{t, x),Il'^{t,y)] = i5'i[5^^\x — y), the creation and annihilation operators appearing 
in f|T4|) should satisfy: 

1- A 



ao(/c),a(j"(/c') 
a||(A?),a+(fc') 



A 

1 - A 
A 



koi2nf6^^\k - k' 
ko{2'Kf6^^\k - k' 



(15) 



For simplicity in the following we will take A = 1/3 so that we use canonically normalized 
operators with positive sign for the temporal component. Notice that this is just the 
opposite situation to the usual Feynman gauge. In fact, A = 1/3 and A = — 1 are 
the only two possible choices with canonical normalizations. As is well-known |19j . 
in order to recover Maxwell's theory, we need to eliminate the negative norm states 
by defining the corresponding restricted Hilbert space. Following the standard Gupta- 
Bleuler formalism, the physical states |0) will be those annihilated by the combination 
ao + ay, that is: (ao(/i;) + a||(fc))|0) = 0. In Minkowski space-time, only transverse degrees 
of freedom contribute to the expectation value of the energy density in the physical states 
and {(j)\Too\(j)) > since the contributions from longitudinal and temporal modes cancel 
each other. Thus, as expected, the theory is free from ghosts. Notice that in Minkowski 
space-time, we also get {(j)\d^A'^\(j)) = 0. 

Now we can proceed to the quantization in the inflationary epoch. In order to 
present the calculational method explicitly, we assume an exact de Sitter phase. The 
general quasi-de Sitter results will be given below. Thus in de-Sitter: a = —l/{Hri) and 
H = —l/rj with rj < 0. The classical solutions of the corresponding equations are a bit 
more complicated, although it is still possible to obtain analytic expressions: 



A 



Ok 



Cikrje' 



kr] 



-{I + ikri)e-'^'^ 



k'^rfe'^'^Ei{2ikri) 



A\\k = iCi{l+ikr])e-'^'^ -iC2 



3 

-e 
2 



■ikr) 



+ (l-ikri)e'^'^Ei(2ikr]) 



(16) 



where Ei[x) = e~^^/tdt is the exponential integral function. Note that the mode 
Ci can be gauged away by means of a residual gauge transformation. The sub-Hubble 
limit {\kri\ ^ 1) of these solutions reads: 



A^ 

A, 



■Ok 



\\k 



-Cikrj — iC2) e 



-ikr} 



(17) 



The choice of adiabatic vacuum [23] is made by matching these solutions with those 
obtained in the Minkowski case ([T3j) (up to sub- leading terms in the sub-Hubble limit). 
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To do this, we choose the modes Ci and C2 in the following way: 
Ci ^ - — (ao + a|j) 

On the other hand, on super-Hubble scales {\kri\ <^ 1) we have: 
Aok = ^C2{kr])-^ 

Ai\k = tCi - - 7 - H'^tkv)] (19) 



2 



OA: 



with 7 the Euler's constant. We see that ^ofe cx a, which means that Aq^ = A 
is (almost) constant during inflation, once the mode leaves the horizon. In the case of 
quasi-de Sitter slow-roll inflation, the Hubble parameter reads H = — 1/((1 — e)?7), where 
the slow-roll parameter is defined as e = l/{167rG){V' /V)"^ <^ 1, with V the inflaton 
potential. Following the same steps as before, we obtain the power spectrum for Aq on 
super-Hubble scales: 



k^ ,, . ,9, H 



2 

2\ _ 



k 



aH, 



riAn 



(20) 



27r^ " ' ' W 

which is almost scale-invariant (as in the scalar field case) since for the electromagnetic 
spectral index we obtain = — 4e. In a similar way it is possible to obtain the 
primordial power spectrum of longitudinal modes on super-Hubble scales: 

(21) 



VAAk) ^ 



IGvr^e^ 



k 



aHj 



(22) 



If we now compare the power spectra for the conformal fields Aq and A\\ we find that: 

which is negligible on super-Hubble scales, and allows us to safely ignore the longitudinal 
modes on such scales after inflation. 

Notice that since e > 0, VAoik) is a red-tilted spectrum which means that the 
contribution to (Aq) from long wavelenghts dominates over small scales. In particular, 
provided inflation lasted for a sufficiently large number of e-folds, this allows to 
decompose the fluctuations field at any given time into a large homogeneous contribution 
(with scales k < TC) and a small inhomogeneous perturbation (fc > H), and therefore 
we can use standard perturbation theory around the homogeneous background. Thus, 
for the homogeneous part we get: 

iAl),.„.= (23) 

where k^,<HQ, N = Ntot — Nq and k^m = e~^ifo is set by the Hubble horizon at the 
beginning of inflation [21]. Here Ntot is the total number of e-folds of inflation which 
should not be confused with Nq which is the number of e-folds since the time when the 
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k f] 

Figure 1. Evolution of Aok and on super-Hubble scales in the radiation 
era. Continuous (dotted) blue lines correspond to Aok for infinite (vanishing) 
conductivity (no difference in the plot). Continuous (dotted) red lines 
correspond to for infinite (vanishing) conductivity. 

scale Hq^ left the horizon. Typical values for Nq are around 50, whereas generically 
there is no upper limit to Ntot- Thus as expected, up to tilt corrections, Hj sets the 
scale for the field dispersion. 

Once the fluctuations generated during inflation enter the radiation dominated era, 
it would be in principle possible that the high conductivity of the universe could modify 
their evolution, spoiling the growing behavior of super-Hubble models found in ([8]). 
However since conductivity does not affect the temporal equation, the Aq modes are 
not modified on super-Hubble scales. In Fig.l we show the evolution of super-Hubble 
temporal and longitudinal modes, both for vanishing and infinite conductivity. We see 
that the evolution exactly corresponds to Aok oc rf as expected from ([8]) in the radiation 
era, even in the infinite conductivity case. We see that A\\k is sub-dominant compared to 
Aok) until the modes re-enter the Hubble radius for \kr]\ ~ 1. The result is not sensitive 
to the change of initial conditions. Thus, the only effect of the high conductivity is the 
damping of the electric field, (which is consistent with the strict neutrality of the plasma). 
In particular, this implies ^jj^^. = ikAok, which corresponds to the field evolution shown 
in Fig. 1. Let us emphasize that the vanishing of the electric field does not imply the 
vanishing of the temporal component ^ofc- 

In fact, it is straightforward to show that the value of V^A^^ giving rise to the 
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effective cosmological constant is not affected by the presence of conductivity. Indeed, 
Maxwell's equations in the presence of conserved currents read: 

V^,F^"' + XVVuA'' = J" (24) 

Taking the four-divergence of the equation we get: 

□ (V^A'^) = (25) 

where we have used current conservation V^J^ = 0. Thus, we see that the field 
V f_i,A^ evolves as a free scalar field, and it is therefore constant on super-Hubble scales, 
independently of the presence of external currents. 



4. Gauge invariance on cosmological scales 

A remarkable consequence of the covariant quantization formalism in the context of 
inflationary cosmology (which had not been considered previously) is the breaking of 
gauge invariance on cosmological scales. Indeed, in this formalism, the classical energy- 
momentum tensor depends on the gauge-fixing term. However, in Minkowski space- 
time, when one takes the expectation value of this object in a physical state (that 
belonging to the restricted Hilbert space) the gauge dependence disappears because 
the contributions from the temporal and longitudinal degrees of freedom cancel each 
other. This is so just because in Minkowski space-time, in the restricted Hilbert space, 
the amplitudes of temporal and longitudinal mode solutions are exactly the same (see 
[T9]). However when considering the quantization in an expanding universe important 
differences arise. At short distances, i.e. for sub-Hubble modes, it is easy to see from 
f|T7|) that the same cancellation takes place, as it should be, and the theory is exactly 
the same as in Minkowski space-time. Therefore the energy density does not depend 
on the gauge-fixing term. Nevertheless, when the modes become super-Hubble, it can 
be seen from (fT9l) . that the amplitude of the temporal modes grows in time faster than 
that of the longitudinal ones. This spoils the mentioned cancellation and a net energy 
density results from the A term. 

Notice that in the covariant formalism the four polarizations are always present. 
In Minkowski space-time (or for sub-Hubble modes) only two of them contribute to 
the energy density, however on cosmological scales also the temporal one can have 
observational consequences. As Aq is also a propagating degree of freedom, the gauge- 
fixing term can be seen in the Gupta-Bleuler formalism as a kinetic term for it, and 
therefore the coefficient A can be fixed by the standard normalization of the creation 
and annihilation operators. This effect could not be studied in other (non-covariant) 
formalisms, such as Coulomb gauge quantization, since only transverse polarizations 
would be present in that case. 

To summarize, the presence of the background cosmological electromagnetic field 
breaks U{1)em symmetry on large scales while preserving local (small-scales) invariance. 
This is analogous to the situation with Lorentz symmetry, where the presence of matter 
or radiation in the Universe breaks global Lorentz invariance, but respecting local 
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transformations. In other words, the presence of a non-vanishing cosmological constant 
could be signalling the breakdown of gauge invariance on cosmological scales. Let 
us emphasize that this effect is a consequence of the quantization of electromagnetic 
theory in the covariant formalism and, as discussed above, it does not modify any of 
the physical predictions of Maxwell's theory for laboratory experiments or astrophysical 
observations. As a matter of fact, the electromagnetic interaction has not been tested 
on distance scales larger than 1.3 AU [25] . 

5. Perturbations and consistency 

Despite the fact that the background evolution in the present case is the same as 
in ACDM, the evolution of metric perturbations could be different, thus offering an 
observational way of discriminating between the two models. With this purpose, we have 
calculated the evolution of metric, matter density and electromagnetic perturbations in 
the longitudinal gauge with (^fyfoo = 2a^$, 5gQi = a^Si, 6gij = a'^{2'^6ij—hij), 6 = Spm/pm 
and taking Afj, = A^°'^{r]) + 6Afj,, where as commented before the main contribution to 
A'lf"^{ri) comes from the temporal component. The propagation speeds of scalar, vector 
and tensor perturbations are found to be real and equal to the speed of light, so that the 
theory is classically stable. We have also checked that the theory does not contain ghosts 
and it is therefore stable at the quantum level. On the other hand, using the explicit 
expressions in [26] for the vector-tensor theory of gravity corresponding to the action 
in ([1]), it is possible to see that all the parametrized post-Newtonian (PPN) parameters 
agree with those of General Relativity, i.e. the theory is compatible with all the local 
gravity constraints for any value of the homogeneous background electromagnetic field. 

In Figs 2-3 we plot the evolution of scalar perturbations, satisfying $fc = in 
this theory, and matter density contrast 6k, in both, vanishing and infinite conductivity 
limits. We find that the only relevant deviations with respect to ACDM appear on large 
scales k Hq and that they depend on the primordial spectrum of electromagnetic 
fluctuations. However, for N\nAQ\> 12, the results on the CMB temperature power 
spectrum and evolution of density perturbation are compatible with observations. 
Taking Mj ~ Trh ~ 10^ GeV in this case, with Tj^j the reheating temperature, we 
obtain the correct value of the dark energy density today. In addition, the different 
evolution of $fc with respect to the ACDM model gives rise to a possible discriminating 
contribution to the late-time integrated Sachs- Wolfe effect [27] . 

The presence of large scale electric fields generated by inhomogeneities in the Aq 
field opens also the possibility for the generation of large scale currents which in turn 
could contribute to the presence of magnetic fields with large coherence scales. This 
could shed light on the problem of explaining the origin of cosmological magnetic fields. 
Work is in progress in this direction. 
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Figure 2. Evolution of with k = 3Hq, corresponding to the maximum 
contribution to the ISW effect. Upper(lower) panel with vanishing(infinite) 
conductivity. Continuous blue line for ACDM, dashed red for A^|nAol 
(upper panel) 18 (lower panel) and dotted green for large A^lnyigl. 




Figure 3. Evolution of the matter density contrast 6k with k = WHq and 
infinite conductivity (vanishing conductivity shows no difference with respect 
to ACDM). Plotted curves correspond to the same models as in the lower panel 
of Fig.2. 
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6. Conclusions and discussion 

We have shown that the present phase of accelerated expansion of the universe can 
be explained by the presence of a cosmological electromagnetic field generated during 
inflation. This result not only offers a solution to the problem of establishing the true 
nature of dark energy, but also explains the value of the cosmological constant without 
resorting to new physics. In this scenario the fact that matter and dark energy densities 
coincide today is just a consequence of inflation taking place at the electroweak scale. 
Such a relatively low inflation scale implies also that no cosmological gravity wave 
background is expected to be measurable in future CMB polarization observations. 

Notice also that any vector-tensor theory (not necessarily electromagnetism) whose 
low-energy effective action is given by ([1]) and in which the vector field only interacts 
gravitationally with the rest of particles would provide a natural model for dark energy. 
In fact all the previous models trying to account for the cosmic acceleration are plagued 
by classical or quantum instabilities, fine tuning problems or inconsistencies with Solar 
System experiments. However, in this work we present, for the first time, an explanation 
to the cosmic acceleration with none of the aforementioned problems. 
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